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Abstract 

The quasiparticle time blocking approximation (QTBA) is considered as a model for the descrip- 
tion of excitations in open-shell nuclei. The QTBA is an extension of the quasiparticle random 
phase approximation that includes quasiparticle-phonon coupling. In the present version of the 
QTBA, the pairing correlations are included within the framework of the BCS approximation. 
Thus, in this model, the + spurious states appear, which are caused by the breaking of the sym- 
metry related to the particle-number conservation. In this work, the method is described which 
solves the problem of the + spurious states in the QTBA with the help of the projection technique. 
The method is illustrated by calculations of + excitations in 120 Sn nucleus. 
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I. INTRODUCTION 



Description of nuclear excitations within the models based on the mean-field approach 
faces the known difficulty of appearing the so-called spurious or ghost states (see, e. g., 
[1]). The spurious state arises because of the breaking of some symmetry. In particular, the 
breaking of the translation symmetry leads to the appearance of spurious l - states. In the 
fully self-consistent theory, for instance, in the self-consistent random phase approximation 
(RPA) or in the self-consistent quasiparticle RPA (QRPA), the broken symmetry is restored 
and the spurious states disappear [1-3]. The self-consistency means here a fulfillment of 
some constraints imposed on the mean-field operator and the amplitude of the residual 
interaction. However, it is difficult to achieve the full self-consistency in the models which 
go beyond the (Q)RPA framework by taking into account additional correlations. So, the 
problem of the spurious states becomes relevant again, even if it was solved on the (Q)RPA 
level. In this work, the quasiparticle time blocking approximation (QTBA, see Refs. [4, 5]) 
is considered. The QTBA is a model intended for the description of excitations in open- 
shell nuclei. The model is formulated in terms of the Green function method. Within the 
QTBA, the single-particle continuum, the pairing correlations, and the quasiparticle-phonon 
coupling (QPC) are included. In this sense, the QTBA is an extension of the QRPA. Note 
that all three components of the model mentioned above are necessary to describe the 
fragmentation and the width of the states in open-shell nuclei. In the present version of the 
QTBA, the pairing correlations are included within the framework of the Bardeen-Cooper- 
Schrieffer (BCS) approximation. Thus, in this model, the + spurious states appear, which 
are caused by the breaking of the symmetry related to the particle-number conservation. 
In the QRPA, the + spurious state is eliminated by taking into account the so-called 
dynamical pairing effects (particle-particle channel contributions, see Ref. [6]). However, in 
the QTBA, this is not sufficient because the spurious state is fragmented due to its coupling 
to the "two-quasiparticle®phonon" (2q®phonon) configurations. In this work, the method 
is described which solves the problem of the + spurious states in the QTBA with the help of 
a combination of the so-called subtraction procedure (see [4, 5]) and the projection technique 
described below. The scheme proposed is illustrated by the calculations of + excitations in 
120 Sn nucleus. 
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II. SYSTEM OF THE QTBA EQUATIONS AND ITS MODIFICATIONS 



The basic equation, which has to be solved within the Green function method to calculate 
the physical observables related to the nuclear excitations, is the equation for the effective 
(renormalized) response function R eS (u). It has the same form both in the (Q)RPA and in 
the QTBA. In the symbolic notations it reads (in what follows we will use notations and 
definitions of Ref. [5]): 

R cS (u) = A(cu) - A(co) T R cS (u) (1) 

where A(uS) is a correlated propagator and JF is an amplitude of the effective residual in- 
teraction. In the case of the QRPA, A{uj) reduces to the uncorrelated 2q propagator A(uj). 
In the general case including pairing correlations, the amplitude T can be represented as a 
sum of two terms 

jr = jr(ph) + jt(pp) ( 2 ) 

where the amplitude jF( ph ) represents interaction in the particle-hole (ph) channel and T^ vv ^ 
includes contributions of the interaction both in the particle-particle (pp) and in the hole-hole 
(hh) channels (in the following for brevity we will use the unified term pp channel implying 
also the hh-channel contributions). The response function R cS (u) enables one to calculate 
the strength function S(E) which determines the distribution of the transition strength 
caused by an external field V°. These quantities are related by the following formulas: 

S{E) = —hnIL(E + iA), (3) 

7T 

Tl(u) = ~Tr((eVyR cS (u)(eV )), (4) 

where IT(a>) is the nuclear polarizability, E is an excitation energy, A is a smearing parameter, 
and e is an effective charge operator. 

One of the important questions arising in the QRPA and QTBA calculations is the 
question of completeness of the configuration space. The size of the basis in this space has 
an impact practically on all the calculated quantities. In particular, configurations with 
a particle in the continuum are responsible for the formation of the escape widths of the 
resonances. The well-known method to include these configurations on the RPA level is the 
use of the coordinate representation within the Green function formalism (see Ref. [7]). This 
method is used in the present approach as described in Ref. [5]. However, incorporation of the 
pp-channel contributions in the coordinate representation leads to considerable numerical 
difficulties. At the same time, the pp-channel contributions (dynamical pairing effects) 
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are very important in the calculations of + excitations in the open-shell nuclei, primarily 
because of the problem of the + spurious state. For this reason, in Ref. [8] a combined 
method was developed, which is a modification of the so-called (r, A) representation proposed 
in Ref. [9] for the QRPA problem. Within this method only the ph channel is treated 
in the coordinate space; the dynamical pairing effects are included in the discrete basis 
representation. 

Consider the general case of the QTBA. By taking into account the decomposition (2) 
one can rewrite Eq. (1) in the form 

R eS (cu) = A (ms+pp) (u) - A (ICS+PP) (u) .F (ph) R cS (uj) (5) 

where propagator (u) is a solution of the equation 

^(res+pp) ^ = A ^ _ A ( u j jtCpp) ^(rcs+pp) ^ _ (g) 

In the present work the version of the QTBA is used in which the ground state correlations 
caused by the QPC are neglected. In this case the correlated propagator A{uS) is defined by 
the equation 

A(uj) = A(uj) - A{uj) A(uj) (7) 

where A(u) is the uncorrelated QRPA propagator, 

$(o;) = ^ cs \uj) -$( rcs )(0), (8) 

and $( res ) (u;) is a resonant part of the interaction amplitude responsible for the QPC in the 
QTBA (see Refs. [4, 5] for details). Combining Eqs. (6) and (7) leads to the new equation 

for A( rcs+pp )(u;): 

^(rcs+pp) ^ = ^ _ £( u j + ^(pp) j A (rcs+pp) ^ _ ^ 

As a result we find that the pp-channel contributions can be included by modification 
of the equation for the correlated propagator, i. e. by replacing Eq. (7) by Eq. (9). The 
modification is reduced to the additional term jF( pp ) added to the amplitude 

Notice, however, that in practice Eq. (5) for R cS (u) is solved in the coordinate represen- 
tation (to take into account the single-particle continuum), whereas Eq. (9) is solved in the 
restricted discrete basis representation. This fact greatly simplifies the problem as compared 
with the initial Eq. (1) in which both the ph-channel contribution and the pp-channel one 
are included in the coordinate representation. At the same time, the use of the restricted 
discrete basis representation for the pp channel is fully consistent with the BCS approxima- 
tion in which the gap equation is solved in the same restricted basis. In more detail, this 
modification of the QTBA equations is described in Ref. [8]. 
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III. THE METHOD OF ELIMINATING THE 0+ SPURIOUS STATE 



The general scheme described above ensures that the energy of the + spurious (ghost) 
state is equal to zero both in the QRPA and in the QTBA. Indeed, it is not difficult to 
show that Eq. (CI) of Ref. [5] for the transition amplitudes of the spherically symmetric 
nucleus in the QRPA has a non-zero solution at the energy u q = and at the total angular 
momentum and the parity J q = + if the gap equation (A25) is fulfilled. The explicit form 
of this solution is 

P(X ost) = c\ 12) P il)ri , (10) 

where C is an arbitrary constant, 

P(i)v = V P W , P W = \l 2j'i + 1 u (1) v {1) , 77 = ±1 . (11) 

Here and in the following it is supposed that the equations of the model are written in 
the representation of the single-quasiparticle basis functions in the doubled space which 
(according to the notations of Refs. [4, 5]) are labelled by the composite indices 1 = {[1], rr^} 
where [1] = {(1), ^j, (1) = {r l ,n l , and r] 1 = ±1 is the sign of the quasiparticle energy 

E 1 = ^E^. That is, the symbol "(1)" stands for the set of the single-particle quantum 
numbers excepting the projection of the total angular momentum m 1? v 2 ^ is the occupation 
probability, and = ^Jl — v 2 ^. 

The existence of the solution (10) means that the response function R eS (u) as a solution 
of Eq. (5) in the QRPA has a pole at uj — 0. The same is true for the QTBA since at u> = 
the QTBA equation (5) coincides with the QRPA one owing to the subtraction procedure 
determined by Eq. (8). 

However, there still remains the following problem: in the QTBA the ghost state can 
be fragmented due to its coupling to the 2q®phonon configurations, despite the energy of 
the dominant ghost state is equal to zero. It can lead to the spurious states at low energies 
distorting respective strength functions. In particular, these fragmented spurious states will 
produce non-zero response to the particle-number operator which has to be exactly equal 
to zero in a correct theory (as, for instance, in the QRPA including pp channel that was 
proved by Migdal, see [6]). 

To solve this problem, the following method is proposed. Let us recast Eq. (9) in the 
form 

^(res+pp) ^ = ^(pp) ^ _ ^(pp) ^ ^ A (rcs+pp) ^ ? ^ 
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where propagator A^ v \u) is a solution of the equation: 

(to) = A(u) - A{uj) F^p) i(pp) (a;) . (13) 

Consider the case = + assuming diagonal approximation for the amplitude (which 
is consistent with the gap equation in the BCS approximation, see Appendix C of Ref. [8]). 
In this case one can keep only diagonal parts of the matrix functions in Eq. (13). In the 
explicit form we have 

(2)r,> M = ^(l)r,, (2)rf M ~ \l)v, W H ^{^V "%>)rf" A$tf'',(2)rf M ( 14 ) 

(34)r)"v"' 

where A$l {2)tf (u) = A^ vA22)rf {u), 

to - 



and ^7)^22)^ is defined by Eq. (C8) of Ref. [8]. 



The matrix function A^^ ,^(0;) has two poles at u = (first and second order ones) 
and can be represented in the form: 

A^\u) = ^- + — + A^ TC ^u) , (16) 



where function A^^^iuj) is regular at uo — > 0. From the symmetry properties of the matrix 
function A^?^ ( 2 w( a; ) ^ follows that the matrices and are real, symmetric: 

a ( 2 ) = = a^ T , a™ = a«* = a« T , (17) 

and possess properties: 

a W=fi x aWfj x , = -rj'aWfj* , (18) 

where 

^(2W =< W < W ( 19 ) 

Further, by substituting (16) into Eq. (13) and putting the coefficients at the same powers 
of to to be equal to each other we obtain the following equations for the matrices and 

a( 2 ) = -i(0) J^ (pp) a (2) = -a (2) J^ (pp) i(0) , (20) 
( 1 + i(0) ^ (pp) ) a (1) = i(0) fj z a^ , (21) 



a W+aWr) z aW=0, (22) 

where 

*7(1)tj, (2),,' = V 8r,tf 6 (12) ■ ( 23 ) 

Eqs. (13) and (16) uniquely determine the matrices and aP-\ However, strictly- 
speaking, from this it does not follow that a solution of Eqs. (20)-(22) constrained by the 
conditions (17) and (18) is unique. Nevertheless, we will find a particular solution of these 
equations assuming that it is a true answer. Calculations confirm the correctness of this 
assumption. 

As follows from the above analysis (see Eqs. (10) and (11)), the solution of Eqs. (20) can 
be represented in the form: 

a (l) V , (2) V > = a P(2)n> ' ( 24 ) 

where a is a constant and p^ is defined by Eqs. (11). In the symbolic notations we have 

a^ = a \p)(p\. (25) 
Let us search for the matrix a« in the form: 

a (1) = f(|p><£i + ie><pi), (26) 

or in the explicit form: 

a (l)r ? , (2) V > = y ( P(l)r, €(2)rf + £(l)r, P(2)rf ) " ( 27 ) 

In addition, in agreement with Eqs. (11), (18), (19), and (27), we will assume that the 
quantity £^ does not depend on rj, i. e.: 

By substituting (25) and (26) into Eq. (21) we obtain: 

£(i) = ZX n-1 ) (12A2) > ( 29 ) 

(2) 



where 



Here 



V /(2j 1 + l)(2j 2 + l) 

Q (12) = ^(1) 5 (12) + i U (l) ~ V (l) ) i U (2) ~ V (2) ) -7(11,22) • (3°) 

POD 

•^fl2,34) = 

*/ 



R^(r) is the radial part of the single-particle wave function, the quantity jF^(r) determines 
the effective residual interaction in the pp channel (see [5, 8]). 
Substitution of Eqs. (25) and (26) into Eq. (22) yields: 

ao 1 = ( 32 ) 
(i) 

Thus, the matrices and are completely determined. 
Let us introduce the following matrices: 

P=l + a^fi z , = l + r} g aW . (33) 

In the explicit form we have: 

P dH { 2W = Krf 5 (i2) + y ( V P(i)^( 2 ) ff + £(i)P( 2 ) ) , (34) 

p Ih (2 ),' = Kv' 5 (i2) + y ( P(i) ?( 2) + V P (2 ) V ) • (35) 
In the general case of non-diagonal matrices these formulas acquire the form: 

P (12)r], (34)t/ = ^r?,r?' ^(13) ^(24) + ^(12) ^(34) a (l)r ? , (3)r?' ^ (36) 
- P (12)» ) , (34)t/' = <W ^(13) ^(24) + ^(12) ^(34) ^ a (l)^, (3)r?' ' ( 37 ) 

It is easy to show that the following equalities are fulfilled: 

PP = P, ptpt = P t ( 38 ) 

Consequently, P and pt are the projection operators. In addition, we have: 

P«M= = a^pt. (39) 

Eqs. (16) and (39) mean that one can use the operators P and Pt to eliminate the coupling 
of complex (2q®phonon) configurations to the ghost state. From the above analysis it follows 
that the elimination can be achieved with the help of the replacement of the amplitude $(o>) 
in formula (9) by the projected amplitude 5> proj (u;) defined as 

$ proj (c;) = P t $(cj)P. (40) 

In what follows this method will be referred to as the QTBA with + projection. 
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IV. CALCULATIONS OF THE 0+ EXCITATIONS IN THE 120 Sn NUCLEUS 



As an illustration of the method described above consider results of calculations of the + 
excitations in the semi-magic nucleus 120 Sn. Calculational scheme is described in detail in 
Ref. [8]. It is based on the Hartree-Fock and BCS approximations and is fully self-consistent 
on the RPA level. The self-consistent mean field and the effective residual interaction (in- 
cluding the spin-orbital and the Coulomb contributions in both quantities) were derived from 
the Skyrme energy functional. In the present calculations, the T5 Skyrme force parametriza- 
tion [10] was used. The single-particle continuum was included completely on the RPA level. 
The set of the phonons entering 2q®phonon configurations in the QTBA calculations in- 
cluded 29 collective modes with values of the spin L in the interval 2 ^ L ^ 9 and with 
natural parity n = (— 1) L . The smearing parameter A in Eq. (3) is equal to 200 keV in all 
the calculations. 
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FIG. 1: Strength function of the isoscalar giant monopole resonance in 120 Sn calculated within the 
QRPA with (solid line) and without (dashed line) contributions of the pp channel. 
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According to this scheme, the following sorts of calculations were performed: (i) the 
QRPA and the QTBA calculations with and without taking into account contributions of 
the pp channel (they will be referred to as pp+ and pp— calculations, respectively); (ii) 
the QTBA calculations with and without + projection (they will be referred to as proj + 
and proj— calculations, respectively). In Figs. 1 and 2, the results of calculations of the 
isoscalar giant monopole resonance are shown. In Fig. 1, the results obtained within the 
QRPA (pp+) and the QRPA (pp— ) are presented. Spurious + excitations arise at low 
energies in the response calculated without contributions of the pp channel. These spurious 
excitations disappear when the pp-channel contributions are included. The same picture 
is obtained in the QTBA calculations. However, to eliminate spurious excitations in the 
QTBA it is insufficient to include the pp channel. Only the QTBA (pp+, proj+) that 
includes pp-channel contributions together with + projection gives the + response free 




Energy (MeV) 

FIG. 2: Strength function of the isoscalar giant monopole resonance in 120 Sn calculated within the 
QTBA including contributions of the pp channel. The results with and without + projection are 
represented by the solid line and the dashed line, respectively. 
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from the spurious states. In the calculations within QTBA (pp+, proj—) the fragmented 
ghost states remain in the low-energy region, though their strength is strongly suppressed 
as compared with calculations in which the pp-channel contributions are not included (see 
Fig. 2). 

In Fig. 3, the + response to the particle-number operator in 120 Sn is shown. These 
results demonstrate that the QRPA (pp— ) and the QTBA (pp+, proj—) produce non-zero 
response. In the calculations within the QRPA (pp+) and the QTBA (pp+, proj+) this 
response disappears within calculational accuracy as it should be. 



i 8 
> 

CD 



S6 
o 

• rH 
■+-> 

o 

a 

Pa 



CD 
-P 



2 - 







1 


120 Sn, + 


II 
II 


particle— number operator 


II 


A = 200 keV 


II 
II 


HF+BCS 


1 1 

l" 
l v 1 

1 1 
1 1 
1 1 
1 1 
1 1 
1 1 


T5 Skyrme force 


QTBA, pp + , proj - 
QRPA, pp - 


1 1 

1 1 









10 15 20 

Energy (MeV) 



25 



30 



FIG. 3: Strength function of the spurious excitations in 120 Sn caused by the particle-number 
operator. The solid line represents calculation within the QTBA including contributions of the 
pp channel, but without + projection. The dashed line represents calculation within the QRPA 
without contributions of the pp channel. 
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V. CONCLUSIONS 



In this work, the projection method is formulated which solves the problem of the + 
spurious states in the model intended for the description of excitations in open-shell nuclei 
with taking into account the single-particle continuum, the pairing correlations, and the 
quasiparticle-phonon coupling. The efficiency of the method is illustrated by the calculations 
of the + excitations in the 120 Sn nucleus within the self-consistent scheme based on the 
Skyrme-Hartree-Fock approximation. 
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